OKHEP-96-04 
Casimir Energy for a Spherical Cavity in a Dielectric: 
Applications to Sonoluminescence 



Kimball A. Milton* 

Department of Physics and Astronomy, The University of Oklahoma, Norman OK 73019 USA 

Y. Jack Ngt 

Institute of Field Physics, Department of Physics and Astronomy, University of North Carolina, 



cn ■ Chapel Hill, NC 27599 USA 

(N 

(February 1, 2008) 

T 1 I 

> 

■ Abstract 

O 

' In the final few years of his hfe, Julian Schwinger proposed that the "dy- 
namical Casimir effect" might provide the driving force behind the puzzling 

r~| ! phenomenon of sonoluminescence. Motivated by that exciting suggestion, we 

. ^ , have computed the static Casimir energy of a spherical cavity in an otherwise 

X. 

H uniform material. As expected the result is divergent; yet a plausible finite 

answer is extracted, in the leading uniform asymptotic approximation. This 
result agrees with that found using ^-function regularization. Numerically, 
we find far too small an energy to account for the large burst of photons seen 
in sonoluminescence. If the divergent result is retained, it is of the wrong 
sign to drive the effect. Dispersion does not resolve this contradiction. In 
the static approximation, the Fresnel drag term is zero; on the other hand, 
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electrostriction could be comparable to the Casimir term. It is argued that 
this adiabatic approximation to the dynamical Casimir effect should be quite 
accurate. 



I. INTRODUCTION 

In a series of papers in the last three years of his life, Julian Schwinger proposed ^ that 
the dynamical Casimir effect could provide the energy that drives the copious production of 
photons in the puzzling phenomenon of sonoluminescence In fact, however, he guessed 

an approximate (static) formula for the Casimir energy of a spherical bubble in water, based 
on a general, but incomplete, analysis [^. He apparently was unaware that one of us had, 
in the late 1970's, completed the analysis of the Casimir force for a dielectric ball 0]. It is 
our purpose here to carry out the very straightforward calculation for the complementary 
situation, for a cavity in an infinite dielectric medium.[| In fact, we will consider the general 
case of spherical region, of radius a, having permittivity e' and permeability yu', surrounded 
by an infinite medium of permittivity e and permeability /x. 

Of course, this calculation is not directly relevant to sonoluminescence, which is anything 
but static. It is offered as only a preliminary step, but it should give an idea of the orders 
of magnitude of the energies involved. It is a significant improvement over the crude esti- 
mation used in Attempts at dynamical calculations exist H^TO|; but they are subject to 
possibly serious methodological objections, some of which will be discussed below. (Other 



theoretical models to explain sonoluminescence are given in |jTTl.) In fact, we anticipate that 
because the relevant scale of the electromagnetic Casimir effect is in the optical region, with 
characteristic time scale t ~ 10~^^s, and the scale of the bubble collapse is of order r ~ 10"^s 
(more relevant may be the duration of each flash, which is ^ lO^^^s), the adiabatic approx- 
imation of treating the bubble as static for calculating the Casimir energy should be very 



preliminary version of this paper appeared in ||7|. 
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accurate. Sonoluminescence aside, this calculation is of interest for its own sake, as one of 



a relatively few nontrivial Casimir calculations with nonplanar boundaries [p!2|-[T9|]. It rep- 
resents a significant generalization on the calculation of Brevik and Kolbenstvedt |2^, who 
consider the same geometry with fie = fi'e' = 1, a special case, possibly relevant to hadronic 
physics, in which the result is unambiguously finite. It is, as noted above, a straightforward 
generalization of the result in ^j; the most significant technical improvement is that here the 
energy is calculated directly. We also examine Fresnel drag and electrostriction; the latter 
may be numerically significant. 

In the next section we review the Green's dyadic formalism we shall employ, and compute 
the Green's functions in this case for the TE and TM modes. Then, in Section 3, we compute 
the force on the cavity from the discontinuity of the stress tensor. The energy is computed 
similarly in Section 4, and the expected relation between stress and energy is found. Fresnel 
drag, in the static approximation, is considered in Section 5, and electrostriction in Section 
6. Estimates in Section 7 show that the Casimir energy so constructed, even with physically 
required subtractions, and including both interior and exterior contributions, is divergent, 
but that if one supplies a plausible contact term, a finite result (at least in leading approx- 
imation) follows. This finite result agrees with that found using (^-function regularization. 
(Physically, we expect that the divergence is regulated by including dispersion.) Numerical 
estimates of both the divergent and finite terms are given in the conclusion, and comparison 
is made with the calculations of Schwinger and others. A simple estimate is given which 
suggests that any macroscopic electromagnetic phenomenon such as the Casimir effect can- 
not possibly supply the energy required for sonoluminescence. There, and in the Appendix, 
where we discuss the form of the force on the surface due to the fluctuating electric and 
magnetic fields, a comparison with M is made. 
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II. GREEN'S DYADIC FORMULATION 



We follow closely the formulation given in P^,|6|. We start with Maxwell's equations in 
rationalized units, with a polarization source P: (in the following we set c = h = 1) 

d 

-VxE = — B, V-B = 0, (2.1) 
where, for an homogeneous, isotropic, nondispersive medium 

D = eE, B = fxB.. (2.2) 
We define a Green's dyadic T by 

E(r, t)= J {dr') dt' r(r, t; r', t') ■ P(r', t') (2.3) 

and introduce a Fourier transform in time 

r°o duj 



r r, 



t-r\t')= r|^e-(*-')r(r,r», (2.4) 



where in the following the uj argument will be suppressed. Maxwell's equations then become 
(which define $) 

V X r = zcj$, V • $ = 0, 

Iv X $ = -iujeT\ V • r = 0, (2.5) 
/i 

in which V = r + 1/e, where 1 includes a spatial delta function. The two solenoidal Green's 
dyadics given here satisfy the following second-order equations: 

(V^ + cj2e^)r' = --V X (V X 1), (2.6a) 

e 

(V^ + a;2e/i)$ = ia;/iV X 1. (2.6b) 



They can be expanded in vector spherical harmonics 121,22| defined by 



— ; =^'Ylm, (2.7) 

//(/ + 1) 



as follows: 



r'(r,r') = lfi{r,r')XU^) + -^V x r')X^^(^]) 



Im 



*(r,r') = {Ury)y^imm - - V X /,(r,r')X,^(f^) 



When these are substituted in Maxwell's equations ( p.5|) we obtain, first 

9i = ~9u fi = h + --Mr-r')y.Un'), 



and then the second-order equations 



(A + uj^lie)gi{r, v') = lujfi J dil" Xl^n") ■ V" x 1, 

(A + ^Ve)/Kr, r) = / rf^^" XL(^^") ■ V" x (V" x 1) 



e 



where 



2d 1(1 + 1) 
Di = Tr^ + --- 



These equations can be solved in terms of Green's functions satisfying 
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which have the form 



Fi{r, r') 



{Di + uj'efi)Fi{r, r') = --5(r - r') 



ik'ii{k'r^)[hi{k'ry) - v4j;(A;V>)], r,r' < a, 
ikhi{kr^)[ii{kr<) — Bhi{kr^)], r,r' > a, 



where 

k = \u\^/JIe, k' = \uj\\J ^'e' , 
and hi = hf^^ is the spherical Hankel function of the first kind. Specifically, 

/,(r,r')=a;VKr,r')X;^(fi'), 



where Fi and Gi are Green's functions of the form (|2.13|) with the constants A and B 
determined by the boundary conditions given below. Given F/, G/, the fundamental Green's 
dyadic is given by 

r(r,r') = 5:LvKr,r')x,„(^])x;^(^]') 



Irri 



(2.16) 



Now we consider a sphere of radius a centered at the origin, with properties e', /x' in the 
interior and e, n outside. Because of the boundary conditions that 

1. 



E_|_, eEj.^ Br, 



B 



(2.17) 



be continuous at r = a, we find for the constants A and B in the two Green's functions in 



Ap 
Bp 
Ag 
Bg 



i/e/Z ei (x' 


)e'i{x] 


- V 


^F]Iei{x)e'i{x') 






A; 

- V 


^¥JIsi{x)s'i{x') 




)e[{x] 


A; 
- V 


^ei{x)e'i{x') 






Ai 
- V 


7J7si{x)s'i{x') 






Ai 




— — IvQjj Ob — 


= k'a, 


the Riccati-Bessel functions 


(x) = xhi{ 


x), 


Sl 


(x) = xji{x), 



(2.18a) 
(2.18b) 
(2.18c) 
(2.18d) 



(2.19) 



and the denominators 



Ai = ^Jen'si{x')e'i{x) - sj e' ^is'i{x')ei{x) , 
Ai = Je'iJ,si{x')e[{x) - JeiJ,'s[{x')ei{x), 



and have denoted differentiation with respect to the argument by a prime. 
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(2.20) 



III. STRESS ON THE SPHERE 



We can calculate the stress (force per unit area) on the sphere by computing the discon- 
tinuity of the radial-radial component of the stress tensor (see the Appendix) : 



T — Tj.r{a—) — Trr{a+), 



(3.1) 



where 



T 



rr--{[e{El-El)-rlJi{Hl-Hl)]). 



(3.2) 



The vacuum expectation values of the product of field strengths are given directly by the 
Green's dyadics computed in Section 2: 



i(E(r)E(r')) = r(r,r'), 

1 



z(B(r)B(r')) 



V X r(r,r') X V, 



(3.3a) 
(3.3b) 



where here and in the following we ignore 5 functions because we are interested in the limit 
as r' V. It is then rather immediate to find for the stress on the sphere (the limit t' ^ t 
is assumed) 



2ia2 i-oo 27r ^ in 



x|(e'-6) 
+ ()"' - )") 



-a^FKa+, a+) + + --r-r'^ G^r, r') 

fj, \ jji' jiar or I 



r=r'=a+ 



r=r'=a+ 



20^7-00 277 An dx 



(3.4a) 
(3.4b) 



where y = ua, 6 = {t — t')/a, and 



InA^A^^ln (5i(xOe;(x) - 5KxOei(x)) - ^(5^^^ jeK^^) + jeK^^)) + constant. (3.5) 



Here the parameter ^ is 



£.J£ — 1 



^4 + 1 



(3.6) 



This is not yet the answer. We must remove the term which would be present if either 



medium filled all space (the same was done in the case of parallel dielectrics ||23|). The 
corresponding Green's function is 



ik'ji{k'r^)hi{k'ry), r,r' < a 
ikji{kr^)hi{kry), r,r' > a 



(3.7) 



The resulting stress is 



D ZTT 



duj 



-iujT ^ ^ 
1=1 



~ 2/ + 1 



Atx 



ei{x)sl{x)] 



x[s'i{x)e\{x) - ei{x)s'l{x)] \. 



(3.8) 



The final formula for the stress is obtained by subtracting ( p.8|) from ( 3.4b ): 



1 ^^.,.g2^(^iL,„A,A, 

2a4 i_oo 27r ^ ^ ' 



1=1 



At[ dx 



+2x'[s[{x')e[{x') - e,{x')s'l{x')] - 2x[s\{x)e\{x) - e,{x)s'l{x)] 



(3.9) 



where we have now performed a Euclidean rotation, 



y ^ iy, X ^ ix, t = t — t' ^ i(x4 — x'^) [S = {x^ — x'^)/a], 
si{x) si{x) = y^//+i/2(x), ei{x) -> e/(x) = ^y^K/+i/2(x). 



(3.10) 



IV. TOTAL ENERGY 

In a similar way we can directly calculate the Casimir energy of the configuration, starting 
from the energy density 

U = i^l±^. (4.1) 



In terms of the Green's dyadic, the total energy is 
E = J{dr)U 
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2i J 

- / ^^e-*-(*^* )y(2/+l) / r^dr 

X i^2k'mr,r) + Gi{r,r)] + l|-r + G,](r, r')|,,=,| , (4.2b) 

where there is no exphcit appearance of e or fi. (However, the value of k depends on which 
medium we are in.) As in |T^ we can easily show that the total derivative term integrates 
to zero. We are left with 

E = - — e-^-^^(2Z + l) / r^dr2e[Fi{r,r) + Giir,r)]. (4.3) 

2l J-oo 2lT Jo 

However, again we should subtract off that contribution which the formalism would give if 
either medium filled all space. That means we should replace Fi and Gi by 

-ik'AF,Gjiik'r)ji{k'r'), r,r' < a 
Fi,Gi = l (4.4) 

I —ikBF,Ghi{kr)hi{kr'), r,r' > a 

so then ( ^TSf ) says 

E = -^(2/ + l) / ^e-'-U / r'drk'%AF + AG)jf{k'r) 

J~oc ZTT I JO 

+ J\^drk^{BF + BG)h^{kr)y (4.5) 

The radial integrals may be done by using the following indefinite integral for any spherical 
Bessel function ji: 



eTrr(r,r) ^TrV x r(r,r)x V 



(4.2a) 



dxx^if{x) = |[((xj/)')^ - jiixji)' - xji{xji)"]. (4.6) 

But we must remember to add the contribution of the total derivative term in ( [4.2b| ) which 
no longer vanishes when the replacement ( [4.4|) is made. The result is precisely that expected 
from the stress (|3.9| ), 

where the derivative is the naive one, that is, the cutoff 5 has no effect on the derivative. 



It is useful here to make contact with the formahsm introduced by Schwinger 0. In terms 
of an imaginary frequency ( and a parameter w, he derives the following simple formula for 
the energy from the proper-time formalism 

E = / dC / dwTi.G, (4.8) 

2tt Jo Jo 

where the trace refers to space, the Green's function is 

G = (4.9) 

w + H 

and the Hamiltonian appropriate to the two modes is (for a nonmagnetic material) 

I TE : doedo - 
H=\ (4.10) 

[ TM : d^-V ■ (l/e)V 

Consider the TE part (the TM part is similar, but not explicitly considered by 
Schwinger). In terms of Green's function satisfying ( p.l2[ ), we have 



1 roo POO °^ ^oo 

E = — dC rfwV(2/ + l) / drr^Fi{r,r-Xh + w), (4.11) 
2n Jo Jo Jo 

where the third argument of the Green's function reflects the substitution in ( p.l2| ) of cu^e — > 
— C^e — w. We now introduce polar coordinates by writing 

(^e + w = p'^, dCdw = ^2p'^ cose dp de, (4.12) 



and integrate over 6 from to 27r. The result coincides with the first term in (4.3). 



V. FRESNEL DRAG 



As may easily inferred from Pauli's book |24[, the nonrelativistic effect of material motion 
of the dielectric, /3(r), is given by the so-called Fresnel drag term, 

E' = J [dr)^-t^(3 . (D X H) = j (rfr)(e/i - l)/3 ■ (E x H). (5.1) 

To preserve spherical symmetry (of course, this is likely not to be realistic) we consider 
purely radial velocities, 
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f3 = (3r. (5.2) 
Then, what we seek is the asymmetrical structure 

f ■ (E(r) X H(r')) = -r ■ (H(r') x E(r)) = --€,,,f, ■ <l>,,(r', r) 



iUl 



Im 



+ [V X Fi{r',r)XU^')] x XL(^^)}. (5.3) 

This is easily seen to reduce to 

r ■ (E X H) = uj--^rJ2Gi{r',r):XiU^') ■ XL(r]) 
r or , 

- ^l;^/ E ^Kr', r)l^imm ■ XL(^^), (5.4) 

Im 

SO when Vl and VL' are identified, and the angular integral is carried out, we obtain the 

corresponding energy for a slow, adiabatic, radially symmetric motion, 

roo du °° (id 1 (9 

1=1 



E' = (3 rHr{e^i-l) / e'^-^ ^(2/ + 1) -^rG^r', r) - -^r'F^r', r) 

(5.5) 

It is clear, immediately, that if the cutoff r is set equal to zero, this vanishes because the 



integrand is odd in compare to (|4.3|) . Since the sign of r is certainly irrelevant, we 



therefore claim that in this quasistatic approximation Fresnel drag is absent. 



Related is the Abraham value of the field momentum p5| , ^ | 



G = ExH, (5.6) 

which gives then an extra contribution to the force density, 

f' = (e-l)|(ExH). (5.7) 

However, as Brevik noted , this is also zero, because, in the Fourier transform, successive 
action of the time derivative brings down uj and —uj. So the continuing controversy about 
which field momentum to use is without consequence here. This should be already obvious, 
because the energy is well-defined, and we have already seen that the force is related to the 
energy by (p7| ). See also the Appendix. 
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VI. ELECTROSTRICTION 



When a dielectric medium is deformed, there is an additional contribution to the force 
density, that of electrostriction [^,0 , 



1, 
2 



fES = I E'p 



(6.1) 



where p is the density of the medium. This term is without effect for computation of the 
force on the dielectric, because it is a total derivative, yet here, where we are calculating the 
stress on the surface, it can be significant. The simplest model for describing the density 
dependence of the dielectric constant is that given by the Clausius-Mossotti equation, 

e-1 



e + 2 



Kp, 



(6.2) 



where i^T is a constant. Consequently the logarithmic derivative appearing in (|6.1j ) is 

p|4(. + 2). 



(6.3) 



The calculation of the electrostrictive Casimir effect for a dielectric ball is given by Brevik 
26| . We have confirmed his result, and generalized it to the situation at hand. Again, the 



contribution if either medium fills all space has been subtracted. The result for the integrated 
stress on the spherical cavity, after the Euclidean transformation is performed, is 



1 °0 POO 



,iy5 



1=1 

(e'-l)(e' + 2) 



2x' (^'^-^'+1/2 (^') 



+ 



(e-l)(e + 2) /2^2 



TT 



- x{Bf + Bg) /^f+i/2(0 + xBg jKly.iO 



(6.4) 
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VII. ASYMPTOTIC ANALYSIS AND NUMERICAL RESULTS 



The result for the stress (p.9|) is an immediate generahzation of that given in 0, and 
therefore, the asymptotic analysis given there can be applied nearly unchanged. The result 
for the energy is new, and seems not to have been recognized earlier. 

We first remark on the special case i/e/I = i/e'/x'. Then x = x' and the energy reduces 

to 

1 /-oo °o J 

E = -^J ^ dy^''' E(2^ + 1)^^ - ^'(Mn (7.1) 



1=1 



where 



If ,^ = 1 we recover the case of a perfectly conducting spherical shell, treated in for 
which E is finite. In fact ( [7.1D is finite for all ^, and if we use the leading uniform asymptotic 
approximation for the Bessel functions we obtain 



Further analysis of this special case is given by Brevik and Kolbenstvedt . 

In general, using the uniform asymptotic behavior, with x = uz, u = I + 1/2, and, for 
simplicity looking at the large z behavior, we have 



E 1 ^—Yj^^Tdz e--VV^^i. In 

27ra ^Jejl ^ J-oo dz 



2 



16^4 I ey ' ^ 



2^ 



(7.4) 



d 

> I dze'^''"'^^^z 

which exhibits a cubic divergence as 5 — > 0. To be more explicit, let us content ourselves 
with with the case when e — 1, e' — 1 are both small and fi = fi' = 1. Then, the leading u 
term is 

E ^> u- dze"^^'^z—- — 

16™ 2J-00 dz{l + z^y 

(7.5) 



64a 4/ 256a 
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Here, the last arguable step is made plausible by noting that since 5 = r/a the divergent term 
represents a contribution to the surface tension on the bubble, which should be cancelled 
by a suitably chosen counter term (contact term). This argument is given somewhat more 



weight by the discussion in ||2^. In essence, justification is provided there for the use of 



zeta-function regularization, which directly gives the finite part here: 

32™ f^^ 2 64a I 4'' ' ^ > 



because Z^^^o — (2 ~ 1)C(~'5) vanishes at s = 2. 

Alternatively, one could argue that dispersion should be included [p8| -|30||, crudely mod- 
elled by 



If this rendered the expression for the stress finite [we consider the stress, not the energy, 
for it is not necessary to consider the dispersive factor d{ijje{ijj)) / duj there], we could drop 
the cutoff 5 and the sign of the force would be positive: (at last, we set e' = 1) 

^ ^ +T28^ ^ Loo (TTIT OT^^' ^ ^ 

where zq = ujoa/v. As z/ ^ 00, zq —>■ 0, and the integral here approaches ttzq/2, and so 

(60 - 1)^ - (eo - 1)^ 3 ,70, 

if we take as the cutofiQ of the angular momentum sum Vc ~ oj^a. The corresponding energy 
is obtained by integrating — 47ra^jF, 

which is of the form of (|7.5| ) with 1/5 oj^a/A. 



^Inconsistently, for then zq ~ 1. If zq = 1 in ( |7.8| ), however, the same angular momentum cutoff 
gives 5/12 of the value in (^^). 
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It is rather more difficult to extract numerical results from the formula for electrostriction, 
Indeed, Brevik |25[ only considers two special cases, e 1, appropriate to a perfect 



conductor, and e — 1 ^ 1. In fact, in the latter case, corresponding to ( [7.5| ), he was able 
to consider only the single I = 1 term in the sum. This is highly unreliable, as such a 
term may be completely unrepresentative (such as having the wrong sign Because 
this electrostrictive stress presents new and somewhat difficult to understand divergences, 
we will defer its consideration to a later publication, and only remark that it is highly likely 
to contribute a term comparable to the finite Casimir estimate. 



VIII. CONCLUSIONS 

So finally, what can we say about sonoluminescence? To calibrate our remarks, let us 
recall (a simplified version of) the argument of Schwinger [|l|]. On the basis of a provocative 
but incomplete analysis he argued that a bubble (e' = 1) in water (e ~ (4/3)^) possessed a 
positive Casimir energy]^ 

^ 47ra3 f (dk) 1 / 1 \ a^K^ A 1 \ 

Ec / ^^-k 1 F ~ 1 F , 8.1 

where K is a. wavenumber cutoff. Putting in his estimate, a ~ 4 x 10~3cm, ~ 2 x lO^cm"^ 
(in the UV), we find a large Casimir energy, Ec ~ 13 MeV, and something like 3 million 
photons would be hberated if the bubble collapsed. 

What does our full (albeit static) calculation say? If we believe the subtracted result, the 



last form in (|7.5|) , and say that the bubble collapses from an initial radius = 4 x 10~3cm 
to a final radius = 4 x 10~^cm, as suggested by experiments f^, we find that the change 
in the Casimir energy is AE ~ +10~^eV. This is far to small to account for the observed 
emission. 



^Note, for small e — 1, Schwinger's result goes like (e — 1), rather than (e — 1)^, indicating that he 
had not removed the "vacuum" contribution corresponding to ( p. 81 ). This is the essential physical 
reason for the discrepancy between his results and ours. 
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On the other hand, perhaps we should retain the divergent result, and put in reasonable 
cutoffs. If we do so, we have 

E = ^' o?K^ ~ -4 X lO^eV, (8.2) 

perhaps of acceptable magnitude, but of the wrong sign. (The emission occurs during the 
collapse.) The same conclusion follows if one uses dispersion, as ( [7.101 ) shows. 

So we are unable to see how the Casimir effect could possibly supply energy relevant 
to the copious emission of light seen in sonoluminescence. Of course, dynamical effects 
could change this conclusion, but elementary arguments suggest that this is impossible 
unless ultrarelativistic velocities are achieved. For example, consider the Larmor formula, 
appropriate to dipole radiation; it gives the power radiated: 

P-^, (8.3) 

where d is the dipole moment. If our bubble, with N atoms, coherently emits radiation, we 
expect 

d ~ 8.4 

where da is an atomic dipole moment, and r a characteristic collapse time for the bubble. 
Thus the energy emitted during one collapse of a bubble in water is 

2 . , / X2 



cmy / \CT]^ 



(8.5) 



So with a ~ 10~^ cm, r ~ 10~^ s (suggested by experiments 0), and da ~ 10~^e-cm, we get 
an energy of only E ~ 10"^'' eV. This is in spite of the assumption of coherent radiation. 
Note that in (|875|), r would have to be ~ 10^^^ s (which is the upper bound to the observed 
flash duration) to yield a total energy of 10 MeV; this would correspond to a velocity across 
the bubble of 10^ cm/s, well in excess of the speed of sound, thus precluding the presumed 
coherent radiation process. 
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We therefore believe that in Eberlein's calculation ||^ there is an implicit assumption of 
superluminal velocities.^ We note that the short wavelength result of Eberlein, (4.7) of the 
first reference in or (10) of the second, can be cast in the dipole form ( ^.3| ) by integrating 
by parts. Up to factors nearly equal to one, 

I d\ a , , 

a-, (8.6) 



where a{t) is the bubble radius. Because a/c < 1, we find that emission energy of 10 MeV 
requires a time scale of te ~ 10^^^ s. This seems to us an implausibly short scale unless 
remarkable relativistic phenomena are involved. [Incidentally, the magnitude of our cutoff 
estimate, (|8.2|) , also agrees with (|8.3|) if if ~ 1/r. This demonstrates that there is nothing 
classical about the estimate (|8.5|)]. 

The only plausible origin of such short time scales lies in the formation of a shock. Such 



a model was sketched by Greenspan and Nadim in |[TT[]. But this picture has nothing to do 
with quantum radiation. 

Recent experimental results have made it even more difficult to accommodate any ex- 
planation based on macroscopic considerations. In particular, Hiller and Putterman |^ 
have found a remarkably strong isotope effect when water (II2O) is replaced by heavy water 
(D2O), where the dielectric properties change by no more than 10%. This, together with the 
already known strong temperature dependence, and strong dependence on gas concentration 
and the gas mixture, may rule out Casimir effect explanations entirely. Yet the subject of 
vacuum energy is sufficiently subtle that surprises could be in store. 



^Indeed, if one follows Eberlein and uses 7 ~ 1 fs (though the experimental value seems to be 
closer to 10 ps) in her model profile, one finds the maximum speed of the bubble surface to exceed 
the speed of light by almost two orders of magnitude. Actually, even with such a small 7, we 
find her result yields an energy output of only 10^'^ MeV, insufficient to explain sonoluminescence. 
Eberlein agrees with these points, and no longer stands by this model pl[ . 
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APPENDIX: DISCUSSION OF FORM OF FORCE ON SURFACE 

There seems, in the literature, to be some confusion about the correct form for the stress 
on a surface due to electromagnetic fields (here, of course, we are interested in vacuum 
expectation values of those fields). The definitive discussion seems to be given in Stratton 



. We also refer here to the manuscript in process by one of us [32|. 
In the text, we computed the force on the surface by considering the discontinuity of the 
stress tensor. 

Tun = \e{El - El) + - Hi), (Al) 

across the surface, where n denotes the direction normal to the surface, and _L directions 
tangential to the surface. This follows directly from a consideration of the interpretation to 
Tnn as the flow of nth component of momentum in the direction n. Because of the boundary 
conditions that 

Ex, D„, Hx, Bn (A2) 
be continuous, the stress on the surface is 



T = T 

•J n 



(A3) 
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in terms of fields on the surface, and where a prime denotes quantities on the "— " side of 
the surface. This is obviously equivalent to the following form for the force density, 

/i/ 



(A4) 



which is just what is obtained from the stationary principle for the energy [32 



The controversy seems to center around the additional "Abraham" term 

f = (e-l)|(ExH). 



(A5) 



(Henceforward we restrict ourselves to nonmagnetic material, /i = 1.) As noted in Section 
5, this makes no contribution to Casimir effect, because the vacuum expectation value is 
stationary. Furthermore, the existence of such a term is dependent upon the (essentially 
arbitrary) split between field and particle momentum. The Minkowski choice for field mo- 
mentum 



G 



M 



D X H 



(A6) 



would not imply this additional force density. The analysis of experimental data given by 
Brevik [^, however, seems to favor the Abraham value. 

If we were calculating a net force on the surface, ( |A5| ) would indeed give a further 
contribution to the force beyond that given by ( |A3| ) . Through use of Maxwell's equations, 
we easily find 

1 



f" 

J n 



e-1 



2e^"^' - \^nE^ + ■ (B5„) + V ■ (EEO 



(A7) 



If n were a fixed direction, the volume integral of this force density would turn into a surface 
integral, and the result given by Eberlein follows. 



F' 



dS 



l]{Bl-Bl) + {e-e')El 



(A8) 



But this result cannot be used to compute the stress. Thus, the formula (C5) given in 
Appendix C of the first paper in is wrong, and, accordingly, so is (3.18) there.0 



^The first derivation there is based on the formula for the force given in the following paragraph, 
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Finally, we note there is yet another formula for the force on a dielctric given in terms 
of polarization charges and currents, 



F = J{dr) 



where 



PpoiE H — jpoi X B 

c 



d 



(A9) 



Ppoi = -V • P, jpoi = — P + cV X M, 



with the polarization and magnetization fields given by 



P = D-E=(1--)D, M = B-H= (/i-l)H. 



(AlO) 



(All) 



Again, it is easy to show that if one is calculating the force in a fixed direction, so one can 
freely integrate by parts, for a nonmagnetic medium, we recover the expected force including 
the Abraham term: 



7712 1 a 

._V. + --(.-l)ExB 



(A12) 



But the integrand in ( [A9p is not interpretable as a force density from which the stress may 
be computed. In effect, it is that interpretation that P| uses. 



while the second is based on an obviously incorrect extrapolation from the vacuum stress tensor, 
which of course gives vanishing stress. 

20 



REFERENCES 



[1] J. Schwinger, Proc. Natl. Acad. Sci. USA 90, 958, 2105, 4505, 7285 (1993); 91, 6473 
(1994). 

[2] H. Frenzel and H. Schultes, Z. Phys. Chem., Abt. B 27, 421 (1934). 

[3] D. F. Gaitan, L. A. Cmm, C. C. Church, and R. A. Roy, J. Acoust. Soc. Am. 91, 3166 
(1992); B. P. Barber and S. J. Putterman, Nature 352, 318 (1991); B. P. Barber, R. 
Hiller, K. Arisaka, H. Fetterman, and S. J. Putterman, J. Acoust. Soc. Am. 91, 3061 
(1992); R. G. Holt, D. F. Gaitan, A. A. Atchley, and J. Holzfuss, Phys. Rev. Lett. 72, 
1376 (1994); R. Hiller, S. J. Putterman, and B. P. Barber, Phys. Rev. Lett. 69, 1182 
(1992); B. P. Barber and S. J. Putterman, Phys. Rev. Lett. 69, 3839 (1992). 

[4] R. A. Hiller and S. J. Putterman, Phys. Rev. Lett. 75, 3549 (1995). 

[5] J. Schwinger, Lett. Math. Phys. 24, 59, 227 (1992); Proc. Natl. Acad. Sci. USA 89, 
4091, 11118 (1992). 

[6] K. A. Milton, Ann. Phys. (N.Y.) 127, 49 (1980). 

[7] K. A. Milton, in the Proceedings of the Third Workshop on Quantum Field Theory 
Under the Influence of External Conditions, Leipzig, September 18-22, 1995, edited by 
M. Bordag (Teubner, Stuttgart, 1966), p. 13. 

[8] E. Sassaroh, Y. N. Srivastava, and A. Widom, Phys. Rev. D 50, 1027 (1994). 

[9] C. Eberlein, Phys. Rev. A 53, 2772 (1996); Phys. Rev. Lett. 76, 3842 (1996). 

[10] A. Chodos, to appear in Orbis Scientiae 1996, Miami Beach, January 25-28, 1996, 
|hep-ph/9604368 . 

[11] E. B. Flint and K. S. Suslik, J. Am. Chem. Soc. Ill, 6987 (1989); C. C. Wu and P. 
H. Roberts, Phys. Rev. Lett. 70, 3424 (1993); H. P. Greenspan and A. Nadim, Phys. 
Fluids A 5, 1065 (1993); M. P. Brenner, R. R. Rosales, S. Hilgenfeldt, and D. Lohse, 

21 



MIT preprint (April 1996). 
[12] T. H. Boyer, Phys. Rev. 174, 1764 (1968); B. Davis, J. Math. Phys. 13, 1324 (1972). 
[13] R. Balian and B. Duplantier, Ann. Phys. (N.Y.) 112, 165 (1978). 
[14] K. A. Milton, L. L. DeRaad, Jr., and J. Schwinger, Ann. Phys. (N.Y.) 115, 388 (1978). 
[15] L. L. DeRaad, Jr. and K. A. Milton, Ann. Phys. (N.Y.) 136, 229 (1981). 
[16] I. Brevik and G. H. Nyland, Ann. Phys. (N.Y.) 230, 321 (1994). 
[17] K. A. Milton, Ann. Phys. (N.Y.) 150, 432 (1983). 
[18] K. A. Milton and Y. J. Ng, Phys. Rev. D 46, 842 (1992). 
[19] C. M. Bender and K. A. Milton, Phys. Rev. D 50, 6547 (1994). 

[20] I. Brevik and H. Kolbenstvedt, Phys. Rev. D 25, 1731 (1982); Phys. Rev. D 26, 1490 
(1983) (E); Ann. Phys. (N.Y.) 143, 179 (1982); 149, 237 (1983). 

[21] J. D. Jackson, Classical Electrodynamics, 2nd ed. (Wiley, New York, 1975), p. 746 ff. 

[22] J. A. Stratton, Electromagnetic Theory (McGraw-Hill, New York, 1941), Chapter 7. 

[23] J. Schwinger, L. L. DeRaad, Jr, and K. A. Milton, Ann. Phys. (N.Y.) 115, 1 (1978). 

[24] W. Pauli, Theory of Relativity (Pergamon, Oxford, 1958), §33. 

[25] I. Brevik, Phys. Rep. 52, 133 (1979). 

[26] I. Brevik, Ann. Phys. (N.Y.), 138, 36 (1970). 

[27] K. A. Milton, Phys. Rev. D 27, 439 (1983). 

[28] I. Brevik and G. EinevoU, Phys. Rev. D 37, 2977 (1988); I. Brevik and I. Clausen, Phys. 
Rev. D 39, 603 (1989). 

[29] P. Candelas, Ann. Phys. (N.Y.) 143, 241 (1982); 167, 257 (1986). 



22 



[30] I. Brevik, H. Skurdal, and R. Sollie, J. Phys. A 27, 6853 (1994). This paper also repro- 
duces the resuhs of 0, but because of a factor of 2 error, fails to note the cancellation 
of the leading e — 1 term. We thank I. Brevik for a discussion on this point. 

[31] C. Eberlein, private communication. 

[32] J. Schwinger, L. L. DeRaad, Jr., K. A. Milton, and W.-y. Tsai, Classical Electrodynamics 
(preliminary version of manuscript in circulation since 1979), now being prepared for 
publication by Addison- Wesley. 



23 



